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一 种 改进 的 多 通道 盲 信号 解 相关 算法 * 
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摘 要 : 针对 多 通道 解 相关 算法 在 混合 信号 协 方差 矩阵 的 最 大 特征 值 较 大 或 者 最 大 与 最 小 特征 值 的 比值 很 
大 (大 于 108 的 病态 混合 ) 时 收敛 速度 变 慢 且 收敛 误差 增 大 的 不 足 ， 本 文 引入 归 一 化 方法 ， 即 在 迭 
代 过 程 中 对 解 混 信号 进行 归 一 化 ， 从 而 限制 其 协 方差 矩阵 最 大 特征 值 的 取 值 范围 ， 且 降低 最 大 与 
最 小 特征 值 的 比值 ,数值 仿真 表明 ， 改 进 后 的 算法 在 降低 达 代 误差 和 加 强 稳定 性 方面 有 了 明显 的 
改善 ， 收 敛 速度 也 得 到 了 可 观 的 提升 ， 改 进 算法 也 具有 更 广 的 适应 范围 . 
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1 引言 


近年 来 ， 育 信和 号 分 离 在 通信 信号 处 理 、 生 物 医学 信号 处 理 、 语 音信 和 号 处 理 等 领域 得 到 了 广 
泛 的 应 用 H3， 一 些 常用 于 言 分 离 的 自 适应 算法 ， 例 如 ， 独 立成 分 分 析 ， 要 求 对 混合 信号 进行 
预 白化 ( 球 化 )， 一 个 零 均 值 随机 矢量 称 为 白 的 ， 是 指 其 协 方差 阵 是 一 个 单位 矩阵 . 在 白化 时 ， 
信号 矢量 z(%k) 用 下 列 变换 进行 预 处 理 


y(k) = Wa(k), (1) 


目的 是 选择 矩阵 WW， 使 得 信号 y = Wz 的 协 方差 矩阵 {fy(X)y(k)”} 变 成 单位 阵 四 ， 白 化 预 处 
理 在 自 适 应 滤波 以 及 多 层 神 经 网 络 等 系统 中 也 常 被 用 作 预 处 理 过 程 来 提高 自 适应 性 能 ， 因 为 它 
能 简化 计算 复杂 度 ， 在 独立 成 分 分 析 中 ， 它 能 减少 近 一 半 未 知 参数 的 计算 四， 

通常 情况 下 ， 多 通道 解 相关 算法 的 预 白化 过 程 稳定 且 快 速 高 效 所 ?然而 ， 当 混合 信号 协 方 
差 矩阵 的 最 大 特征 值 较 大 以 及 最 大 与 最 小 特征 值 比值 很 大 (大 于 109 的 病态 混合 ) 时 ， 会 出 现 收 
敛 性 能 恶化 的 现象 . 针对 这 一 问题 ， 我 们 在 多 通道 解 相关 算法 中 引入 了 对 每 一 步 解 混 信 号 的 归 
一 化 ， 从 而 控制 了 最 大 特征 值 的 取 值 范围 ， 降 低 了 最 大 与 最 小 特征 值 的 比值 ， 这 样 改进 后 的 算 
法 在 降低 迭代 误差 和 加 强 稳 定性 方面 有 很 大 的 提高 ， 同 时 收敛 速度 提高 了 至 少 一 个 数量 级 . 

本 文 第 2 部 分 简要 介绍 了 全 局 和 局 部 自 适应 两 种 多 通道 解 相关 算法 ; 第 3、4 部 分 分 析 了 算 
法 的 性 能 ， 并 针对 原 算 法 的 不 足 提 出 了 一 种 改进 方法 ， 另 外 还 给 出 了 改进 算法 的 步 长 系数 的 经 
验 公式 ; 最 后 一 部 分 是 改进 前 后 算法 的 数值 仿真 对 比 . 


收 稿 日 期 : 2008-04-10. 作者 简介 : 徐 仔 波 (1984 年 7 月 生 )， 男 ， 硕 士 . 研究 方向 : 育 信 号 处 理 及 其 在 地 震 信号 处 
理 中 的 应 用 . 
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2 独立 成 分 分 析 与 多 通道 解 相 关 算 法 


独立 成 分 分 析 在 盲 信 号 分 离 中 占据 着 重要 的 地 位 ， 其 混合 模型 为 z = 4s 十 vw， 其 中 z 为 混 
合 信号 ，4 为 未 知 混合 矩阵 ，s = [s , sn] 了 为 相互 独立 的 源 信号 ，w 为 观测 噪声 。 其 分 离 模 
型 为 y = Wa, RPW ADAE. 一般 还 假设 s 的 各 个 分 量 的 均值 为 零 、 方 差 为 1 且 4 为 列 
满 秩 矩 阵 等 . 它 利 用 一 个 神经 网 络 全 来 对 混合 信号 > 进行 处 理 ， 通 过 无 监督 的 自 适应 学 习 来 
优化 一 个 对 照 函数 (Contrast function)， 并 不 断 调整 网 络 权 值 WW， 使 得 W 和 混合 矩阵 4 的 乘 
积 逐 渐 收 敛 至 一 个 广义 排列 矩阵 互 ， 此 时 网 络 输出 y 即 为 源 信 号 波形 的 估计 值 回 .在 上 述 过 程 
中 ， 白 化 预 处 理 过 程 显得 尤为 重要 ， 而 常用 的 预 白化 处 理 算法 就 是 多 通道 解 相关 算法 . 

典型 的 多 通道 解 相关 算法 有 Silva 和 Almeidal5:6] 所 提出 的 全 局 自 适应 算法 





W(k+1)=W(k) +n(k)(T — y(k)y(K)7)W (k), (2) 

以 及 Cichocki, Kasprzak 和 Amari 提 出 的 局 部 自 适应 算法 [8 
W(k 4-1) = W(k) + n(k) [I — y(5)y(5)7]. (3) 
两 种 算法 都 涉及 一 个 问题 ， 如 何 选 取 W(0) 和 wm(k)， 从 而 使 算法 有 最 佳 的 收敛 性 能 ?Douglas 和 
Cichocki 给 出 了 上 述 两 种 算法 的 误差 演变 式 以 及 步 长 系数 的 选取 方法 ， 使 得 系数 矩阵 W(k) 能 

够 快速 、 稳 定 、 稳 健 地 收敛 四. 
对 于 全 局 自 适应 算法 外 有 

Ki(k+1)= [(1 — n(k)Ai(k))? — n(K)? A0] Ac). (4) 


上 式 是 对 (2) 式 的 误差 演变 分 析 , HEH ilk) = Ae) — 1, Alk) 为 Ry) 的 第 i 个 特征 值 , Ry, (k) 
B x SL D, Sc EA [0] 中 式 (19), nk) 为 步 长 系数 . 由 (4) Ñ, Douglas 和 Cichocki 得 出 了 保证 每 
^r Alk) 稳定 地 收敛 到 1 的 、n(k) 应 满足 的 充分 条 件 


0 < n(k) < Mmax(k), (5) 


若 0 Amax(k) <3+2 2, 
m=] Jm # O< Xmox(k) «32/2 " 


Xo Ë 和 max(k) > 3 十 2V2， 
此 即 文献 g] 中 的 (22) 与 (23) 两 式 ， 式 中 as (R) 为 最 大 步 长 系数 ， 和 wmax(k) 为 Ruy (k) KBAR 
Ec 自 适应 算法 (3)， 见 文献 [9] 的 式 (28)， 有 
walk +1) = [1 — n() Xa t XQ )] Au), (7) 
上 式 是 对 算法 (3) 的 误差 演变 分 析 ， 其 中 A 和 wi(k) = Au (E) -1/ Ai XE ON W (k) RÆ, Au (K) 


为 BE{W(k)} 的 第 i 个 特征 值 ， 和 i 为 Ras 的 第 i 个 特征 值 ，Ai(k) = 和 iX2,;(k). 利用 式 (7); Doug- 
las I Cichocki F H T RUE Aw a(k) 收敛 到 1/V 和 的 、n(k) 应 满足 的 稳定 性 条 件 


0<n(k) < (8) 


2 
VAmax(1 + VAmax(k)) 
其 中 Amax 和 Aca (k) 分 别 为 Rss M Ryg (c) 的 最 大 特征 值 . 
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3 ”算法 性 能 分 析 


当 上 述 Amax Ñ Amax(k) 很 大 时 ， 由 式 (6) 和 式 (8) TA, BKRS N) 取 值 变 得 很 小 ， 此 
时 收敛 速度 变 慢 ， 如 图 1(a)、 图 2(a) 所 示 . 特别 地 ， 当 混合 信号 协 方差 矩阵 的 最 大 特征 
f Dass 与 最 小 特征 值 和 win 之 比 (定义 A = Amax/Xmin) 很 大 (大 于 108 的 病态 混合 ) 时 ， 若 初 
fü W(0) = I, BE Ry, (5) 初始 变化 不 大 ， 即 RSS) = Rss， 所 以 误差 演化 式 (4) 和 式 (7) 可 以 
近似 为 





(+) > | 位 - 2:9 *- 389,0) = ailh), 


入 max 


Kailh + D) e [1 - 8/502 09]5,, i(k) = 83,9), 


Amax 


(9) 


显然 ， 对 于 那些 小 特征 值 和 A; 之 Amn 的 误差 传 符 项 a 和 6B (a 和 BB 由 上 式 确定 ) 将 接近 1， 此 
时 ，Ryy(k) 将 停留 在 距离 单位 矩阵 了 的 一 定 误差 处 而 不 收敛 至 I， 如 图 3 所 示 . 













x 最 大 特征 值 =1215.2 
—— 最 大 特征 值 =247.0 
-最 大 特征 值 =22.4 


一 -~ 最 大 特征 值 =1215.2 
co 最 大 特征 值 =247.0 
c 最 大 特征 值 =22.4 





10. 
一 一 Tum E 最 大 特征 值 =3459.0 
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(a) 改进 前 收敛 图 (b) 改进 后 收敛 图 


Bi: 改进 前 后 的 全 局 自 适 应 算法 应 用 于 不 同 最 大 特征 值 的 混合 信号 收敛 图 
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(a) 改进 前 收敛 图 (b) 改进 后 收敛 图 


图 2: 改进 前 后 的 局 部 自 适应 算法 应 用 于 不 同 最 大 特征 值 的 混合 信号 收敛 图 
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(a) 全 局 自 适应 学 习 算法 (b) 局 部 自 适应 学 习 算 法 


图 3: ”改进 前 后 全 局 和 局 部 自 适 应 算法 应 用 于 人; = 1.022 x 107, Amax = 35076 的 混合 信号 收敛 图 


4 ”多 通道 解 相关 改进 算法 
我 们 对 于 现 有 的 解 相关 算法 提出 如 下 改进 : 在 每 一 步 解 相 关 迭 代 之 前 ， 将 解 混 信号 y 进行 
妇 一 化 
z(k) = G(k)y(k), (10) 
其 中 P : 
G - dag([ S MD y (mop) ^). (11) 
t=1 t=1 


即 解 混 信号 y 在 时 域 上 的 归 一 化 矩阵 ，diag(.) RRR () 中 元 素 为 对 角 元 素 的 对 角 和 矩阵 ，7 为 
信号 时 域 长 度 ，n 为 信号 维 数 ， 显 然 ， 归 一 化 对 于 y 中 的 大 功率 (大 特征 值 ) 分 量 有 压制 作用 ， 
而 对 于 小 功率 (小 特征 值 ) 分 量 有 放大 效果 ， 即 归 一 化 前 后 Az(k) € Ay(k)， 加 快 了 收敛 速度 . 
此 时 ，z 的 每 一 行 元 素 的 平方 和 为 1，Rss(k) = E(z(k)z(k)7] = z(k)z(k)” 是 一 个 对 称 正定 矩 
阵 ， 主 对 角 元 素 为 1，RRss(i,i) = 1 > |Rsz(i,7)|， 其 中 i,j = 1,…,n 且 i 关 7. 所 以 ， 对 于 具 
有 大 的 和 max(k) 的 思 jy(k)， 上 述 变 换 将 最 大 特征 值 降 到 一 个 合适 的 水 平 可 以 证 明 ， 当 z 的 维 
BUS ni, R.Q(k) 的 最 大 特征 值 和 max(k) 满足 0 < Az max(k) < n. 这 可 以 解决 收敛 速度 随 最 
大 特征 值 增 大 而 减 慢 的 问题 ， 且 由 于 最 大 特征 值 和 ,max(k) € (0,n)， 根 据 式 (6)、 式 (8)， 步 长 
系数 nm(k) 的 取 值 范围 (0, +00) 将 变 成 (+é, 十 co)， 而 的 取 值 跟 信号 维 数 有 关 . TE, KRA 
上 述 自 适应 学 习 算法 得 到 如 下 改进 算法 


W(k+1)= W(k) + n()[I — GG)y(E)u()* G(R)] W (k), (12) 
W(k +1) = W(k) + n(k)[I — G(k)y(k)u(k)" G(K)], (13) 


其 中 y(k) = W(k)z(k), GWR (11) 所 示 . 
通过 大 量 的 数值 仿真 实验 ， 我 们 给 出 改进 后 的 全 局 和 局 部 自 适应 算法 的 最 佳 步 长 系 
数 mort 的 经 验 公 式 
(n +1) (n/A; ) 全 局 算法 ， 
Tlopt = _1 、 (14) 
(n 十 1)?(2nVAs) , 局 部 算法 ， 


其 中 n 为 信号 维 数 ，A; 为 协 方差 矩阵 Ros 的 最 大 特征 值 与 最 小 特征 值 之 比 . 
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5 ”数值 仿真 


1) ”病态 混合 以 及 不 同 最 大 特征 值 的 混合 信号 的 解 相 关 过 程 

在 Matlab 环 境 下 ， 选 取 源 信和 号 s = [si, 82, 83, S4, 85) H s2 = 1， 取 混合 矩阵 4 为 随机 矩 
阵 进行 混合 的 同时 利用 调幅 参数 a; 改变 各 源 信 号 幅 值 ， 从 而 得 到 不 同 最 大 特征 值 和 wsx 和 
大 As 值 的 混合 信号 T (病态 混合 )， 即 t= A|[a151, 2252, a383, 0484, asss]. 

图 1、 图 2 中 四 条 收敛 曲线 对 应 混合 信号 zx 的 最 大 特征 值 和 ,sx 依次 为 3459.0, 1215.2, 247.0, 
22.4, A, 依次 为 2169,1710, 598, 559. 根据 式 (5), (8) 以 及 式 (14)， 图 1(a)、 图 1(b) 和 图 2(a)、 
图 2(b) 中 最 大 特征 值 和 max 从 大 到 小 对 应 学 习 步 长 系数 n(k) 依次 取 为 0.00056, 0.0016, 0.0081, 
0.061; 0.105, 0.11, 0.22, 0.19; 0.00028, 0.0008, 0.0042, 0.05; 0.06, 0.07, 0.12, 0.10. 由 图 1、 图 2 
可 知 ， 改 进 算法 步 长 系数 n(k) 并 不 随 最 大 特征 值 和 wax 变化 而 急剧 变化 ， 而 是 在 小 范围 内 波 
动 ， 收 敛 稳定 且 快 速 ， 全 局 自 适应 算法 的 收敛 步 数 随 最 大 特征 值 在 100-2500 步 内 变化 ， 最 大 特 
征 值 越 大 收敛 越 慢 ; 改进 后 全 局 自 适应 算法 的 收敛 步 数 大 都 在 20 步 左右 . 局 部 自 适应 算法 的 
收敛 步 数 随 最 大 特征 值 在 250-2500 步 内 变化 ， 最 大 特征 值 越 大 收敛 越 慢 ， 改 进 后 局 部 自 适应 算 
法 的 收敛 步 数 大 都 在 30-50 步 左右 ， 图 3 显示 了 改进 后 的 算法 在 最 大 特征 值 与 最 小 特征 值 的 比 
值 较 大 情况 下 的 收敛 性 能 . 

根据 式 (5), (8) 和 式 (14) 估计 最 佳 学 习 步 长 系数 ， 图 3(a)、 图 3(b) 中 改进 前 后 的 全 局 和 局 
部 自 适应 算法 学 习 步 长 系数 依次 取 为 0.000056, 0.0011 和 0.000025, 0.0007. 显然 ， 改 进 后 的 算 
法 不 会 出 现 收敛 误差 大 的 问题 且 收 敛 速度 加 快 . 全 局 算法 的 收敛 步 数 达到 30000 时 ， 依 然 存 在 
较 大 误差 ， 且 很 难 继续 收敛 至 单位 矩阵 ， 而 改进 后 全 局 算法 大 约 在 3500 步 左右 收敛 至 理想 结 
果 ; 如 图 3 (b) 所 示 ， 局 部 算法 与 全 局 算法 类 似 ， 改 进 的 局 部 算法 大 约 在 5000 步 左右 收敛 至 理 
想 结果 . 

2) 运算 时 间 的 分 析 与 统计 

局 部 自 适 应 算法 : 原 算法 中 存在 2 次 矩阵 相 乘 ; 在 改进 算法 中 ， 存 在 3 次 矩阵 相 乘 ， 这 
里 G(k) 的 时 间 计 算 复杂 度 可 以 忽略 不 计 . 全 局 自 适 应 学 习 算法 : 原 算法 中 存在 3 次 矩阵 
相 乘 ， 而 改进 算法 中 存在 4 次 矩阵 相 乘 . 因此 ， 对 于 局 部 自 适应 算法 ， 改 进 后 每 一 步 的 时 
间 计 算 复 杂 度 大 概 提 高 了 50%， 而 全 局 算法 大 概 仅 提高 了 33 双 ， 然 而 ， 友 代步 数 的 减少 量 
在 10 ~ 103 量 级 ， 单 步 的 时 间 复 杂 度 对 于 整个 算法 的 运算 速度 的 影响 就 可 以 名 咯 不 计 了 。 

计算 环境 如 下 : 计算 机 型 号 为 HP Compaq de 7800 Convertible Minitower; 处 理 器 型 
55 Jj Intel(R) Core(TM)2 Quad CPU Q6600 92.40GHz; A # (RAM) 为 4.00GB， 得 到 了 如 
表 1 所 示 的 统计 结果 . 虽然 在 单 步 运算 时 间 上 改进 算法 更 长 ， 但 整体 耗 时 要 比 原 算 法 少 ， 实 际 
运算 时 间 与 上 述 的 分 析 基 本 一 致 ( 表 中 全 局 算法 的 平均 每 步 耗 时 稍 小 于 局 部 算法 的 平均 每 步 耗 
时 ， 可 能 与 计算 机 的 稳定 性 有 关 ). 特别 地 ， 对 于 某 些 大 特征 值 混合 信号 的 解 相关 过 程 ， 改 进 
算法 运算 时 间 远 远 比 原 算 法 少 ， 例 如 对 于 Amax 为 3459.0 时， 如 表 1 TR. 





表 1: 不 同 解 相关 算法 的 运算 时 间 统 计 








特征 值 3459.00 1215.20 247.00 22.40 | 平均 步 时 

局 部 算法 2745 步 Tcal 步 /0.17s | 55435 /0.063s | 52 37 /0.0061s | 0.000098s 
全 局 算法 2633 步 /0.22s | 1182 步 /0.12s | 248 步 /0.029s | 56 步 /0.0067s i 0.000092 s 
改进 局 部 算法 | 53 步 /0.014s | 43 步 /0.011s | 27 步 /0.0073s | 28 步 /0.0075s | 0.00027s 
改进 全 局 算法 | 22 步 /0.0060s | 22 步 /0.0060s | 14 步 /0.0040 E 35 /0.0038 d 0.00027 s 
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6 结论 


本 文 在 全 局 和 局 部 两 种 解 相关 算法 中 引入 了 对 解 混 信 号 的 归 一 化 ， 跟 制 了 解 混 信和 号 协 方差 
和 矩阵 的 特征 值 (控制 其 协 方差 矩阵 的 最 大 特征 值 取 值 范围 且 降 低 了 最 大 与 最 小 特征 值 之 比值 )， 
使 算法 在 降低 误差 和 加 快 收 敛 速度 方面 得 到 了 提高 ,同时 ， 改 进 后 的 算法 具有 相对 稳定 的 收敛 
性 能 ， 能 适应 于 更 广 的 应 用 范围 . 





参考 文献 : 


1 乔 永 凤 ， 马 建 芬 . 独立 分 量 分 析 技术 在 语音 信号 处 理 中 的 应 用 [可 . 电脑 开发 与 应 用 ，2007, 20(4): 22-23 
Qiao Y F, Ma J F. Application of independent component analysis technology in speech signal processing[J]. 
Computer Development & Applications, 2007, 20(4): 22-23 

[2] Hyvärinen A, Oja E. Independent component analysis: algorithms and applications[J]. Neural Networks, 
2000, 13(4-5): 411-430 

[3] Vigario R. Independent component analysis for identification of artifacts in magnetoencephalographic 
recordings[J]. Advances in Neural Information Processing Systems, 1998, 10: 229-235 

[4] Cichocki A, Amari S I. Adaptive Blind Signal and Image Processing[M]. LTD: John Wiley & Sons, 2002 

[5] Silva F M, Almeida L B. A distributed solution for data orthonormalization[C]// International Conference 
of Artificial Neural Networks, Espoo, Finland, 1991, 1: 943-948 

[6] Silva F M, Almeida L B. A distributed decorrelation algorithm[J]. Neural Networks, 1991: 145-163 

[7] Cichocki A, Kasprzak W, Amari S I. Multi-layer neural networks with a local adaptive learning rule for 
blind separation of source signals(C]// International Symposium of Nonlinear Theory Application, Las 
Vegas, NV, 1995: 61-65 

[8] Comon P. Independent component analysis-a new concept?[J]. Signal Processing, 1994, 36: 287-314 

[9] Douglas S C, Cichocki A. Neural networks for blind decorrelation of signals|J]. IEEE Transactions on Signal 
Processing, 1997, 45(10): 2829-2842 


An Improved Algorithm for Blind Multi-channel Signal Decorrelation 
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Abstract: As the maximal eigenvalue of the mixture's covariance matrix becomes very large or the 
ratio of the maximal eigenvalue to its minimal eigenvalue is very large (which is larger than 10° as 
ill-conditioned mixture), the current blind multi-channel signal decorrelation algorithms tend to be 
slow in convergence and large in decorrelation error. In this paper, we propose a new normalization 
algorithm which restricts the range of the covariance matrix's maximal eigenvalue and depresses the 
ratio of the maximal eigenvalue to its minimal eigenvalue. The new algorithm can greatly improve the 
convergence speed, decorrelation error, and computational stability of the blind multi-channel signal 
decorrelation. The simulations verify these advantages of the proposed algorithm. 

Keywords: decorrelation; eigenvalue; ill-conditioned mixture; normalization; blind signal separation; 
independent component analysis 
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